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ABSTRACT 

We study the behavior of the three-point correlation function ^ of dark matter 
and mock galaxies, concentrating on the effects of redshift-space distortions and the 
determination of galaxy bias parameters in current redshift galaxy surveys. On large 
scales, redshift space distortions tend to wash out slightly the configuration dependence 
of the reduced 3-point function Q3 ^ C/C^- On smaller scales (< 10/i~^Mpc), Q3 
develops a characteristic U-shape anisotropy between elongated and open triangles 
due to the effects of velocity dispersion. We show that this shape is quite universal, 
very weakly dependent on scale, initial spectral index or cosmological parameters and 
should be detectable in current galaxy surveys even if affected by shot-noise or galaxy 
bias. We present a detailed method for obtaining constraints on galaxy bias parameters 
from measurements of Q3 in current galaxy redshift surveys, based on the eigenmode 
analysis similar to the one developed for the bispectrum. We show that our method 
recovers the bias parameters introduced into mock galaxies by a HOD prescription 
and is also able to handle potential systematics in the case when a smaller number 
than ideal of mock catalogs is used to estimate the covariance matrix. We find that 
current redshift surveys (e.g. SDSS or 2dFGRS) are just about large enough to get 
interesting new constraints on bias. 



1 INTRODUCTION 

The galaxy three-point function provides a valuable statis- 
tical tool to investigate the nature of the relationship be- 
tween galaxies and dark matter and to probe the statistical 
properties of primordial fluctuations through constraints on 
primordial non-Gaussianity (see Bernardeau et al. 2002 for 
a review). 

Measurements of the three-point function and other 
higher-order statistics in galaxy catalogs have a rich his- 
tory (Peebles & Groth 1975, Fry & Peebles 1978, Baumgart 
& Fry 1991, Gaztafiaga 1992, Bouchet et al. 1993, Fry & 
Gaztaiiaga 1994). In the past decade, three-point statistics 
have confirmed the basic picture of gravitational instabil- 
ity from Gaussian initial conditions (Frieman & Gaztafiaga 
1994, Jing & Borner 1998, Frieman & Gaztaiiaga 1999, 
Feldman et al. 2001). The impact of such measurements 
on theoretical models has been important but the system- 
atic and statistical uncertainties in the data have been non- 
negligible, primarily because the connection between these 
statistics and theoretical predictions is best done on large 
scales, where the physics is best understood, and the sur- 
veys previously available were not large enough. 

In recent years, with the completion of large redshift 
surveys such as 2dFGRS (CoUess et al. 2001) and SDSS 
(York et al. 2000) the measurement of higher-order statis- 



tics promises to provide tight constraints at the few percent 
level (Colombi et al. 1998, Szapudi et al. 1999, Matarrese 
et al 1997, Scoccimarro et al. 2004, Sefusatti & Scoccimarro 
2004) that will help significantly in pinning down the cos- 
mological parameters and issues in galaxy formation that 
have an impact on galaxy clustering. First measurements of 
three-point statistics in these redshift surveys are given in 
Verde et al. (2002), Jing & Borner 2004, Kayo et al. (2004) 
and Wang etal (2004). 

In this paper we study the three-point function with 
particular emphasis on redshift-space distortions and galaxy 
bias. There is a puzzling result regarding previous investiga- 
tions of the 3-point function in redshift (configuration) space 
(Suto & Matsubara 1994, Matsubara & Suto 1994, Jing & 
Borner 1998, Jing & Borner 2004). They all find very little 
dependence of the 3-point function on triangle configuration 
at small scales (both in simulations and observations). This 
result is puzzling because one of the first things that brings 
attention to our eyes is the presence of the "fingers of god" 
in redshift-space distributions, produced by the radial com- 
ponent of galaxy peculiar velocities inside clusters. These 
rather long and thin structures should leave a clear imprint 
in the 3-point correlation function, with higher amplitudes 
for collapsed triangles. Why has this feature not been ob- 
served so far? 
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We also study the three-point function as a probe of 
galaxy bias in the framework of current surveys, and provide 
a detailed framework for computing reliable error estimates 
on bias parameters that include the covariance matrix due 
to non-lincaritios and survey geometry. Wc use a battery 
of N-body simulations and mock catalogs extracted from 
them with different biasing prescriptions. Our presentation 
follows closely the eigenmode analysis method developed for 
the bispectrum (Scoccimarro 2000). 

A first paper of this series (Barriga & Gaztanaga 2002, 
paper I) presented a comparison of the predictions for the 
two and three-point correlation functions of density fluctu- 
ations, ^ and in gravitational perturbation theory (PT) 
against large Cold Dark Matter (CDM) simulations. Hero 
we extend these results into the non-linear regime and fo- 
cus on the effects of redshift distortions and the extraction 
of galaxy bias parameters in galaxy surveys. This paper is 
organized as follows. In section 2 we review the basic def- 
initions involving the three-point function. In section 3 we 
present the simulations and mock catalogs, whereas in sec- 
tion 4 we study the effects of redshift distortions from large 
to small scales. In section 5 we study in detail how to recover 
bias parameters from current galaxy catalogs. 



2 THE THREE-POINT FUNCTION 
2.1 Definitions 

The two and three-point correlation functions are defined, 
respectively, as 



ari2) = (<5(ri)<5(r2)) 
C(n2,r23,ri3) = {S{ri)5{r2)S{r3)} 



(1) 
(2) 



where 6{r) — p(r)/p— 1 is the local density fluctuation about 
the mean p = (p) , and the expectation value is take over dif- 
ferent realizations of the model or physical process. In prac- 
tice, the expectation value is over different spatial positions 
in our Universe, which are assumed to be a fair sample of 
possible realizations (see Peebles 1980). It is convenient to 
define a Qs parameter as follows (Groth & Peebles 1977) 



Qs = 



C{ri2,r23,ri3) 



(3) 



CH{ri2,r23,ri3) 
Ch = ^(ri2)^(r23)+e(ri2)C(m)+^(r23)^(ri3), (4) 

where we have introduced a definition for the "hierarchi- 
cal" three-point function (h- The Q3 parameter was thought 
to be roughly constant as a function of triangle shape and 
scale (Peebles 1980), a result that is usually referred to as 
the hierarchical scaling. As we will see below, accurate mea- 
surements/predictions show that Q3 is not quite constant in 
any regime of clustering, although the variations of Q3 with 
scale are small compared to the corresponding changes in ^ 
or f . 



2.2 Triangle parametrizations: a vs. v 

There are two parametrizations that are commonly used in 
the literature to describe the three parameters that define 
a triangle, each of which has advantages and disadvantages. 



Since we are studying the monopole of the three-point func- 
tion (averaged over all orientations of the triangle with re- 
spect to the observer), this depends on three variables. One 
obvious choice is the magnitude of the three sides of the 
triangle, described by fi2, r23 and r3i = fi2 + r23. This is 
reasonable since the three variables are on an equal foot- 
ing (being three lengths) and if the measurement is done 
with some typical resolution due to the three-point function 
algorithm, it affects all of them in the same way. 

The two most common parametrizations, however, in- 
volve a length scale, a dimensionless number (ratio of 
lengths) , and some measure of triangle shape (another ratio 
of lengths or an angle) . In this case the three variables may 
be affected by resolution (e.g. binning) in different ways. 
This can be important for example for the dependence on 
triangle shape, which is what concerns us here regarding the 
signature of fingers of God at small scales. 

In perturbation theory, it has been customary to use 
the two sides of the triangle ri2 and r23 (which are typically 
comparable in side) and the angle a between them: 

/ N ri2 r23 fr-s 
cos{a) = . (5) 

rrz r23 

with < a < 180 deg. For small scales it has been popular 
to use ri2 and r23, but with later given by the ratio: 



u = r23/ri2, 

and the third parameter given by: 



r31 - ^23 

ri2 



(6) 



(7) 



with the restriction ri2 < r'23 < ris. The value of v is related 
to a by: 



V = \/ V? — 2 u cos a-\-\ — u 



(8) 



In principle — 1 < w < 1, but the restriction ri2 < r23 < ri3 
only allows < w < 1. The minimum value oi v = corre- 
sponds to cos a — l/{2u). Note how in this parametrization 
the small values of a are only recovered for large values of 
u. For example, for u = 1 we are restricted to 60 < a < 180 
deg. This can have implications for the detection of sharp 
dependencies on triangle shape, as we shall see in section 4. 

Elongated or "collapsed configurations" are those with 
a ~ or a ~ 180 deg; or in terms of v, those with v close 
to unity. Configurations of triangles with a ~ 90 are called 
here "perpendicular" configurations. In terms of v they cor- 
respond to small values of v, but note that there is no one- 
to-one correspondence. Small values of v could also arise for 
triangles with r3i = r23 even if ri2 is not small. 

We use the term "strong configuration dependence" 
when there is a significance difference between the collapsed 
and the perpendicular configurations. By "weak configura- 
tion dependence" we mean that Qs is hierarchical (ie con- 
stant as a function of a). 



2.3 3-point Function Estimator and Algorithm 

We follow the fast algorithm described in some detail in 
paper I. This algorithm allows a fast calculation of two and 
three-point function for millions of points. The first step is 
to discretize the simulation box into Lsize^ cubic cells. We 
assign each particle to a node of this new latticed box using 
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Table 1. Main parameters describing the simulations. Box size Lsize is in h Mpc. For VLS and HV the shape parameter F stands for 
an eflfective value corresponding to the lineaj: power spectrum with a transfer function for Clm = 0.26, f2(, = 0.04 and Ha = 0.7 = h. 



the nearest grid point particle assignment. We preealculate 
the Ust of relative neighbors to any given node in the lattice. 
To compute now the two-point and three-point correlation 
functions we use: 

ari2) = (9) 

C(n2,r23,r-i3) = — (10) 

where i extends over all nodes in the lattice, j over the list of 
precalculated neighbors (see paper I) that are at a distance 
ri2 from i and k over the neighbors at distance riz from j 
and ri3 from i. 

Besides the tests presented in Barriga & Gaztafiaga 
(2002) we have performed a series of additional tests de- 
signed to quantify how robust our new findings are in red- 
shift space and on smaller scales. We have checked that very 
similar results arc obtained for different pixel sizes. Pixels 
can be as large as ri2/3 on a side with negligible effects in 
the resulting Qz(q). This of course imposes a constraint on 
how thick we can take the binning in a and works as long 
as Qa has a smooth shape as a function of a. As we will 
show below (ic Fig. 5) the pixel scale or the binning play an 
important role on small scales. 

Dilution of point density by factors of 1/10 or 1/20 does 
not bias the results, indicating that the above estimator is 
robust. Using realistic mock catalogues (eg. with a 2dFGRS 
or SDSS mask) instead of cubical boxes gives very similar re- 
sults. In addition, the LCDM240 and LCDM400 simulations 
(see table 1 below) give values of Qa in excellent agreement 
at small scales. 



& Gaztafiaga 1996 for more details). The particle mass is 
M ~ 2.3 X lO^^M0//i. 

A second set, which we just call ACDM240 consists 
on a single realization of the ACDM model normalized to 
(Ts = 0.9 in a box L = 239.5ft" ^ Mpc with IS! = 256^ parti- 
cles. This simulation was run with a parallel version of the 
AP3M code by the VIRGO consortium. The particle mass 
is M ~ 6.9 X 10^°Msu„//i. Thus, while a single ACDM400 
has ~ 4 times more volume (20 times more volume if we 
use all 5 realizations) the ACDM240 model has ~ 33 times 
higher particle density (and therefore 33 times higher mass 
resolution). This will allow us to explore how Qa is affected 
by resolution and sampling effects. The difference in normal- 
ization [a% = 1.0 compared to as = 0.9) is quite small and 
has very little impact on Qa. 

We also study throe further different models using sim- 
ulations from the VIRGO consortium with identical volume 
and number of particles as ACDM240. The different mod- 
els arc listed in Table 1. The comparison with the SCDM 
allows us to study the dependence of Qa on the shape pa- 
rameter r, while the rest of the models explore sensitivity 
to cosmological parameters and normalization. 

Finally, we consider two other simulations from the 
VIRGO consortium, the GIF (higher resolution), VLS 
(larger volume) and HV (even larger volume but with lower 
resolution) simulations, also shown in Table 1. These in par- 
ticular were used to build the galaxy mock catalogs by us- 
ing the Halo Occupation Distribution (HOD) formalism (Ma 
and Fry 2000, Peacock and Smith 2000, Scljak 2000, Scocci- 
marro et al. 2001, Berlind et al 2003, see Cooray and Sheth 
2002 for a review), see Section 5.1 for details. 



3 SIMULATIONS 

We use several sets of simulations centered around the 
Lambda Cold Dark Matter model (ACDM from now on) 
with Qa = 0.7, Qm = 0.3, h = 0.7. The first set, which 
we call ACDM400, consist of 5 independent realizations 
of ACDM model normalized to as = 1 in a box L — 
378ft~^Mpc with N = 126^ particles. These simulations 
were run with the AP3M code of Efstathiou et al. (1985) to 
study volume-averaged higher-order correlations (see Baugh 



4 Qa IN REDSHIFT SPACE 

Redshift distortions (Jackson 1972, Sargent and Turner 
1977, Peebles 1980, Kaiser 1987) duo to peculiar velocities 
along the line of sight create an anisotropic galaxy distribu- 
tion with a preferred direction, that joining the observer to 
the each galaxy. Therefore different survey geometries will 
be affected differently by redshift distortions. Here we focus 
on the long-distance observer approximation where wc place 
the observer in one of the axes at infinity. This allows a gen- 
eral study of redshift distortions, independent of the survey 
geometry. Nevertheless we do not find strong differences for 
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Figure 1. Particle positions in real (loft) and redshift space (right) in a 1 Mpc/h thick projected slice from to the LCDM240 simulation 
(see table 1). 
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Figure 2. Q3(ri2, r23, «) as a function a, the angle between ri2 and ^23, in ACDM 400 simulations. Real space measurements are shown 
as dashed lines, while dotted lines show tree-level (real-space) perturbation theory predictions. Symbols with errorbars correspond to 
redshift space measurements. Each panel shows the results for triangles of different scales. 
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more realistic situations, including that of a redshift slice 
(where one of the angular distances, typically declination, is 
restricted to a few degrees on the sky). 

Figure 1 shows a 1 Mpc/h slice in the ACDM240 case. 
One can see that: 

a) fingers of God are long and sharp needles of about 
l^~^Mpc wide and 10/!-~^Mpc long. This should make a 
very noticeable imprint on Qa for elongated triangles (i.e. 
the "collapsed" configurations in §2.2) across those scales. 

b) filamentary structure is apparent on large scales in 
both cases, although in redshift space it seems to be more 
washed away. This should show as weaker configuration de- 
pendence (or less contrast between collapsed and prependic- 
ular triangles) for Qs in redshift space than in real space at 
these large scales. 

c) on smaller (non-linear) scales structures are rounder, 
and seem sharper in real space. Thus, weaker configuration 
dependence and larger values of Qa are expected on small 
scales in real space. 

Let's check how these visual impressions are reflected 
in the actual measurements of Q3. 

Figure 2 compares the estimations of Q3 in the 
ACDM400 case for triangles with |fi2| = |r23|/2 = 
1.5, 3, 6, 9, 12, 18/i~^Mpc in real space (dashed lines) and 
redshift space (symbols with errorbars). The symbols with 
error-bars show the mean and variance in 5 realizations of 
the same ACDM simulations, while the dotted lines show 
the corresponding leading order (tree-level) PT prediction 
in real space. 

On large scales Q3 in redshift space has only slightly less 
configuration (or shape) dependence than in real space, as 
expected from perturbation theory (Hivon et al. 1995) and 
in agreement with bispectrum results (Scoccimarro et al. 
1999), indicating that part of the angular asymmetries are 
washed out by the peculiar velocities, but this effect is small. 
Note that on the largest scales 18 — 36/i~^Mpc, the values 
of Q3 both in real (or redshift space) seem to approach the 
PT predictions, within the large errorbars. This provides an 
important observational test for the current picture of struc- 
ture formation, independent of cosmological parameters for 
a fixed power spectrum shape (see Bernardeau et al. 2002) . 

On scales smaller than < 10/i~^Mpc (where ^ > 1), the 
redshift space results are almost independent of the overall 
scale of the triangle (despite the large changes in the am- 
plitude of (,), while the real space values increase sharply as 
a function of scale (at least up to ~ l^~^Mpc). Note how 
Q3 for collapsed configurations q ~ 0, 180 deg in redshift 
space rises sharply, while for the perpendicular configura- 
tions (q ~ 90) the values of Q3 are quite fiat and scale 
independent. This characteristic anisotropy is what we call 
the U-shape. This U-shape does not appear in real space 
or in PT theory, as it is basically the results of small-scale 
redshift distortions, the effect of fingers of God. A similar ef- 
fect is seen in the small-scale redshift-space bispectrum (see 
Figs. 3 and 4 in Scoccimarro et al. 1999). 

In Fig. 3 we compare the values of Q3 in redshift space 
for different models on small scales. The U-shape anisotropy 
is present in all models regardless of differences in the cos- 
mological parameters. Note also the small errorbars and the 
good agreement between the two ACDM models with very 
different resolution (and slightly different erg normalization). 
These results extend to smaller scales: we find very similar 
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a 

Figure 3. Q3(ri2, ^23, a) in redshift space as a function a, the 
angle between ri2 and r23, for simulations in Table 1. Symbols 
with 3-sigma errorbars (from the scatter in 5 realizations) corre- 
spond to ACDM400. Continuous, short dashed, long dashed and 
dotted lines correspond to OCDM, ACDM240, tCDM and SCDM 
respectively. 



results (with no evidence for any change with scale) in the 
range ri2 — 0.1 — 3 Mpc/h. All cases in Fig. 3 correspond 
to ri2 = r23/2 = l.OMpc/h, except for the ACDM400 (sym- 
bols with errorbars) which has ri2 — t'23/2 — 1.5Mpc/h. 
Figure 4 shows the comparison between ACDM and SCDM 
on different scales. As shown in the figure there is a small 
global shift between the two models on all scales. On the 
largest scales, the difference become smaller (in contrast to 
what happens in real space, see paper I). 

The results of Fig. 3 show that the U-shape anisotropy is 
expected in a wide range of models, and as we will see below 
(see §5.1) it is also expected in biased galaxy distributions. 
However, it has not been seen previously in N-body simu- 
lations or observations. Early measurements in simulations 
(Suto & Matsubara 1994, Matsubara & Suto 1994) claimed 
that the hierarchical ansatz {Q3 constant independent of 
configuration and scale) is valid in redshift space at small 
scales. More recent studies (Jing & Borner 1998) see a very 
mild effect, but nothing comparable to the amplitude seen 
in Fig. 3. In observations, similar results are seen in LCRS 
(Jing & Borner 1998), 2dFGRS (Jing & Borner 2004) and 
SDSS (Kayo et al. 2004). We believe that these discrepancies 
with our results are due to these authors' choice of binning, 
as we discuss now. For example, Matsubara & Suto (1994) 
bin distances which vary by 30%, Jing & Borner (1998,2004) 
have only 5 bins in the v variable, comparable to Kayo et al 
(2004), who use 5 bins in a from to 180 deg. 

In order to investigate possible explanations for this sit- 
uation, we explore the effect of using different pixel resolu- 
tions and sparse sampling of triangles. We see in Fig. 5 that 
this typically results in a biased-low value of Q3, significantly 
so at collapsed configurations. Because there are far fewer 
triangles in the collapsed configuration, low-resolution trans- 
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Figure 4. Comparison of redshift-space Q3{ri2,r23,ct) for 
SCDM400 (tiiick linos) and ACDM400 (thin lines). 

lates into a reduced U-shape, effectively mixing a small num- 
ber of large amplitudes with large numbers of small ampli- 
tudes and thus down-weighting the resulting Q3 at collapsed 
configurations. A similar effect is obtained if the binning in 
angle is made too wide; the effect is even more severe in the 
w- variable parametrization (see section 2.2) which has been 
traditionally used at small scales. We believe these reasons 
may explain why the U-shape anisotropy has not been seen 
before in the three-point correlation function studies. 



5 GALAXY BIASING 

We now check whether biasing affects significantly the char- 
acteristic U-shape dependence we find in redshift space at 
small scales. We will also study what constraints we can get 
on biasing parameters from the measurements of Q3. 

The value of Q3 depends on galaxy bias, in a form that 
can be easily calculated on large enough scales, where galaxy 
fluctuations So can be modeled as a local (but non-linear) 
function of the corresponding matter fluctuations 5; thus 
expanding this local function for 5 ^ 1 it follows that (Fry 

6 Gaztanaga 1993) 

So ^ F[S] ^ E ^ '5*' (11) 

k 

where k — comes from the requirement that (So) ~ 0. 
It then follows (see Fry & Gaztafiaga 1993, Frieman & 
Gaztariaga 1994) that: 

~ 1 {Qs + C2) (12) 
01 

where C2 = &2/&1, and the ~ sign indicates that this is the 
leading order contribution in the expansion given by Eq. (11) 
above. 
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Figure 5. Effects of resolution and sparse sampling. The dashed 
black line shows the result with highest resolution, which gives 
the higher results for the collapsed configurations (o = 0, 180) . 
The lowest line corresponds to the lowest resolution (cubic cells 
of 1 Mpc//i and sparse sampling of 1 inlO triangles). Intermediate 
lines shows an increase in resolution by a factor of two (with full 
sampling). 

Thus, in general, the linear bias prescription is not accu- 
rate for higher-order moments even when S <^ 1, the reason 
being that nonlinearities generate non-Gaussianities of the 
same order as those of gravitational origin. The linear bias 
term &i can produce distortions in the shape of Q3, while the 
non-linear terms C2 only shifts the curve. Therefore, it is pos- 
sible to use the shape of in observations: when compared 
to the dark matter predictions for Qs one can separate 61 
from 62 in the above relation. This gives an estimate of the 
linear bias 61 which is independent of the overall amplitude 
of clustering (eg. erg). This approach has already been im- 
plemented for the skewness S3 (Gaztanaga 1994, Gaztanaga 
& Frieman 1994), the bispectrum (Frieman & Gaztafiaga 
1994, Fry 1994, Feldman et al. 2001, Verde et al. 2002) or 
the angular 3-point function (Frieman & Gaztafiaga 1999). 

We now test these ideas with mock galaxy catalogues. 

5.1 Mock galaxy samples 

We use three sets of mock catalogues, which we will called 
SdFGRS mocks, HV mocks and HOD mocks. Each set has 
different charactertistics which can be used to test the dif- 
ferent approximations used. 

5.1.1 2dFGRS mocks 

The 2dFGRS mocks consist of 22 independent mock 2dF- 
GRS surveys, extracted from the HV simulations (Evrad et 
al 2002) and described in Norberg et al. (2002). They have 
the same radial and angular selection function as the final 
(public) version of the 2dFGRS and have been convolved 
with the completeness mask of the survey (Cole etal 2005). 
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We extracted 4 volume limited samples shown in Table 2. 
These same samples have been used by Baugh et al. (2004) 
and Croton et al. (2004a,b) to study higher-order moments 
of counts-in-cells and void probabilities. Unless stated oth- 
erwise, we will use as a default the volume limited mock 
with —21 < Alb < —20, which yields the best constraints 
(see Fig. 14). 

5.1.2 HV mocks 

These consist of a set of 100 cubic subvolumes of side 
L = 350 Mpc//i, extracted directly from the HV simula- 
tion (see Table 1). There is no galaxy biasing, dilution or 
realistic mask involved in these mock catalogs. We will use 
them to estimate close to ideal conditions in a dark mat- 
ter distribution. The HV mocks each have 2-3 times more 
effective volume than the 2dFGRS mocks. 

5.1.3 HOD mocks 

We also build galaxy mock catalogs following the results 
from the measurements of the two-point function in the 
SDSS (Zehavi et al. 2004). For the thresholds Mr < -20 and 
Mr < —21, we use the VLS simulation, and a halo occupa- 
tion distribution (HOD) prescription for the mean number 
of galaxies in a halo of mass m > mmin 

/ m \ " 

(iVgai(m)) = 1 + exp[-ms/(m - mmin)], (13) 

\mi / 

and zero otherwise, where the first contribution is that due 
to a central galaxy, the rest being satellite galaxies which 
are taken with a Poisson distributed scatter (Kravtsov et al. 
2003). The power-law index a is fixed to 1. For Mr < — 19 
and Mr < —18 mock galaxies were kindly provided by 
A. Berlind, based on the GIF simulation, with the same 
HOD as in Eq. (13) except for a sharp cutoff at mmin in- 
stead of an exponential cutoff, and the power-law index a is 
allowed to vary. The parameters of the HOD's are given in 
Table 3. 

Given a HOD, the expected large-scale bias parameters 
(in the limit of large scales) are given by, 

bi = — dm n{m) {N„^i{m)) bi{m), (14) 
J 

where n(m) is the halo mass function (which we take to be 
that of Sheth and Tormen 2002), bi(m) are the correspond- 
ing halo bias parameters (Scoccimarro et al. 2001) and the 
galaxy number density is given by 

Wg = y dm n{m) {iVgai(m)). (15) 

Note how the different absolute magnitude limited sam- 
ples have different biasing parameters, contrary to what hap- 
pens in the 2dFGRS mocks, which have a fixed biasing pre- 
scription for all magnitude ranges. The main reason for using 
the HOD mocks is that we know the values of hi , 62 that we 
input in our catalogues, and we can therefore test how well 
bi can be recovered. 

Figure 6 shows how well these values of 61 and 62 fit the 
mean biasing scatter relation 5g = F[S] in Eq. (11) for cubes 
of side I = 30 Mpc/h (see also Fig. 1 in Scoccimarro 2000 
for more cases) . As can be seen in the figure, the agreement 




Figure 6. Scatter plot showing galaxy density constrast 5g ver- 
sus dark matter (5jv/ fluctuations in the Mr < —20 HOD sample 
for cubic cells of side Z = 30 Mpc/h. The smooth long-dashed 
line shows the bias relation with parameters bi and 62 given by 
the HOD from Eq. (14) with parameters in Table 3. The true 
mean bias relation is represented by the central thick line and 
68% scatter around it (parallel thin lines). The short-dashed line 
is the unbiased prediction, bi = 1, 62 = 0. 

is quite good. Similar results are found for other smoothing 
scales: I > 15 Mpc/h. On smaller scales, the scatter plot is 
dominated by shot-noise (the number of galaxies per pixel 
becomes too small). 

The HOD mocks do not have a realistic survey mask, 
as in the 2dFGRS mocks, but instead correspond to a sim- 
ple square box geometry, as in the HV mocks. Given the 
large size of the VLS simulation (about 3 times in volume 
the HV mocks) and its simple geometry, a single mock can 
be used as a prediction for the theoretically expected cor- 
relation functions for each set of bias parameters. The 22 
2dFGRS and 100 HV mocks will be used to estimate errors 
and covariance matrices (see §4). 

Figure 7 shows a slice corresponding to the HOD mock 
galaxies with Mr < —20. Because of the low density we 
have used a thicker slice here than in the dark matter case 
of Fig. 1. Note how the fingers of god are not prominent 
in this case and is hardly possible to detect them visually 
(compare with Fig. 1). We now study whether we can still 
detect in Qs the signature of fingers of god, the U-shape 
anisotropy. 

5.2 Galaxy Qa: Testing the validity of the bias 
expansion 

Figure 8 illustrates how biasing changes the values of Q3 
from the dark matter simulations (symbols with errors) to 
the HOD mock galaxies (continuous lines); all quantities 
here are in redshift space. The errors in the dark matter 
simulation are obtained from the 100 HV mocks. Similar re- 
sults are found for other magnitude ranges. Note that the 
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Table 2. Properties of the combined 2dFGRS SGP and NGP volume limited catalogues (VLCs). Columns 1 and 2 give the faint and 
bright absolute magnitude limits that define the sample. Column 3 gives the median magnitude of the sample, computed using the 
Schechter function parameters quoted by Norberg et al. (2002). Columns 4, 5 and 6 give the number of galaxies, the mean number 
density and the mean inter-galaxy separation for each VLC, respectively. Columns 7 and 8 state the redshift boundaries of each sample 
for the nominal apparent magnitude limits of the survey; columns 9 and 10 give the corresponding comoving distances. Finally, column 
11 gives the combined SGP and NGP volume. All distances are comoving and are calculated assuming standard ACDM. 
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-18.0 
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-18.44 


23290 
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0.014 
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0.0056 


0.021 


0.130 


61.1 


375.6 


7.97 


-20.0 


-21.0 


-20.28 


33997 


1.46 
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0.033 


0.188 


95.1 
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23.3 


-21.0 
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-21.16 


6895 


0.110 


0.0209 
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rig 


bi 


C2 


Mr 
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-18 


1.87 X 10" 


3.74 X 10^2 




0.92 


0.026 


0.97 


-0.21 


Mr 


< 


-19 


3.87 X 10" 


8.8 X 10^2 




1.08 


0.013 


1.07 


-0.08 


Mr 


< 


-20 


9.9 X 10" 


1.71 X 10^3 


6.86 X 10^2 


1 


0.005 


1.13 


-0.11 


Mr 


< 


-21 


5.09 X 10^2 


5.77 X 10^=* 


5.77 X 10" 


1 


0.001 


1.32 


-0.03 



Table 3. Parameters of the mean halo occupation distribution used to populate dark matter halos in simulations with mock galaxies, 
see Eq. (13). All masses are in units of Mq/H. The last three columns correspond to the number density rig (in galaxies per Mpc^/i"^) 
and bias parameters calculated from Eq. (14), with = h2lh\. 
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Figure 8. Q3 in redshift space dark matter simulations (points 
with errorbars, from 100 HV mocks) and respective results for 
Mr < —20 HOD mocks (continuous line). Dashed lines shows the 
HOD bias prediction using the values from Table 3 in Eq. (12). 



effects caused by biasing at the smallest scales are com- 
parable to the variations among different models shown in 
Fig. 3. The characteristic U-shape anisotropy due to fingers 
of god is also present in the galaxy samples on small scales 
ri2 < 6/i~^Mpc, even though it is not so obvious to the eye 
in Fig. 7. 

The HOD prediction based on Eqs. (14) and (12) 
(dashed lines) seems to work very well, even approaching 



the non-linear regime the deviations are not too large. Note 
however that the error bars become very small so these devi- 
ations are statistically significant at small scales, as we shall 
see in more detail below (see eg. bottom panel in Fig. 12). 
On the other hand, it is encouraging that the validity of 
the bias expansion seems to be much broader than perhaps 
naively expected, at least for the galaxy HOD's we have 
studied here. In this section and below we assume we can 
trust the local bias expansion, more precisely Eq. (12) for 
quantities in redshift space, where we always use simulations 
to make predictions for Q3 for dark matter in redshift space. 

To characterize how well we can recover the actual bias 
parameters predicted by a given HOD it is necessary as we 
discuss now to study in detail the covariance between mea- 
surements of Qs for different triangle shapes and scales. 

5.3 Q-eigenmodes 

As we shall see in the next section, it is reasonable to expect 
the hkelihood function to be Gaussian in terms of Q3, at 
least as a first approximation in the case of large surveys, 
which means that one can estimate parameters using a 
approach (likelihood £ oc exp— x^/2), 

= E E ^^^T/^. (16) 

i=i 

with Cij the normalized (to unit variance) covariance ma- 
trix and A, EE [Qt'ii) - QT'^{i)]/cTQ{i), where Qf""^ cor- 
responds to the model (i.e. Eq. (12) with simulations used 
to predict the dark matter Q3), is the measured re- 

duced three-point function of the galaxy population whose 
bias parameters we want to determine, and erg are the un- 
certainties in Qs which are typically obtained from mock 
galaxy catalogs (as the diagonal part of the covariance ma- 
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Figure 7. Galaxy positions in the real (left) and redshift space (right) in a 5 Mpc//i slice corresponding to HOD mock galaxies with 
Mr < —20. The observer is at a; — > — oo, just as in Fig. 1. 



trix), or estimated directly from the observations eg. using 
independent sub-samples or jacknife resampling. The label 
i = I, . . . , Nb denotes different triangle shape configurations. 
In our case, without any loss of generality, we have fixed 
ri3 = 2ri2 and bins refer only to changes in a with Nb = 18 
(i.e. Aq — 5 dcg.). In general one could mix many differ- 
ent shape configurations (i.e. different scales) and consider 
a larger covariance matrix (Scoccimarro 2000). 

The normalized covariance matrix is estimated from the 
mock catalogs as. 



— r 



A': A'] 



(17) 



with Ai = [Qsii) — Q3{i)]/i7Q{i), where Qsii) is the mean 
of the Nm mock independent realizations denoted by the 
index k. Figure 9 shows the estimated covariance matrices 
as a function of the configuration angle a in triangles of scale 
ri3 = 2r'i2 = 24 Mpc/h. As can be seen in this figure, the 
covariance is quite significant and cannot be neglected. The 
structure is very similar between the 2dFGRS (top) and HV 
(bottom) cases. 

Before we invert dj in Eq. (16) , notice that the values of 
Cij are estimated in practice to within a limited resolution, 



ACi, 



(18) 



therefore if Nm is small or if there are degeneracies in dj 
the inversion will be affected by numerical instabilities. In 
order to eliminate this problem, we perform a Singular Value 
Decomposition (SVD) of the matrix, 

Cij = {Uik)^DkiVij, (19) 

where Dtj — X^Sij is a diagonal matrix with the singular 
values on the diagonal, and U and V are orthogonal matri- 
ces that span the range and nuUspace of Cij (see Eisenstein 
& Zaldarriaga 2001, for an application to the power spec- 



trum estimation). By doing the SVD decomposition, we can 
choose the number of modes we wish to include in our 
by effectively setting the corresponding inverses of the small 
singular values to zero. In practice, we work only with the 
subspace of "dominant modes" which satisfy: 



A? > ^/2/N^ 



(20) 



which is the resolution to which we can estimate the co- 
variance matrix elements. In other words, when doing the 
inverse of C in Eq. (16), one inverts Eq. (19) with a restric- 
tion on Dki based on the criterion given by Eq. (20). 

Within the Gaussian assumption one needs to go no 
further than what we have discussed to obtain constraints 
on bias parameters. However, Gaussianity must be chocked 
for each particular survey and scale under study. As shown 
for the bispectrum (Scoccimarro 2000), the non-Gaussian 
case can be simplified by rotating into the basis where the 
Qa's become diagonal, and their non-Gaussian distributions 
can be obtained from mock catalogs and multiplied together 
effectively assuming that they are independent to construct 
the likelihood function for the bias parameters. 

In other words, one can represent each of the mock 
or data Qa vectors (the vector index i = 1, Nb = 18 runs 
over different values of a in this case) in a new basis of Q- 
eigenvectors where the covariance matrix becomes diagonal 
(since U is the rotation matrix for nonzero modes ), 



Q3(i) = ^(7,, 



Q3(j) 



(21) 



The data expressed in the new basis, the Q-eigenmodes Qz, 
is a vector which in general has a lower dimensionality than 
Nb due to the criterion given by Eq. (20). Contrary to Qa, 
the Q-cigcnmodcs are effectively statistically independent 
(this becomes rigorous in the Gaussian case), in the sense 
that their covariance matrix is diagonal. The eigenvalues 
Af are the diagonal elements of their covariance matrix, i.e. 
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Figure 9. Normalized covariance matrix dj as a function of the 
configuration angle a (in degrees) from Q3(o) for fixed side tri- 
angles with ri3 = 2ri2 = 24 Mpc//i. White and black colours 
corresponds to +1 and —1 covariance, respectively. Middle grey 
corresponds to zero values. The top panel shows the covariance 
matrix from 22 2dFGRS mocks, while the bottom panel is ob- 
tained from the 100 HV mocks. 



the variance of a given Q-eigenmode. One can construct a 
"signal to noise" (S/N) for each eigenmode labeled by i, 




The total S/N can be obtained by adding the individual 
modes in quadrature. When we use only the dominant Q- 
eigenmodes, the total S/N estimate corresponds to a lower 
bound, which is not optimal, but avoids potential errors of 
a singular inversion. Figure 10 shows the three most im- 
portant (in terms of signal to noise) Q-eigenvectors for the 
case of the HV mocks at scales ria = 2ri2 = 24 Mpc//i. 
As can be seen in the figure the Q-eigenvectors correspond 
to interesting linear combination of different triangles: the 
first mode (solid) is basically the average Q3 (slightly down- 
weighted for collapsed configurations, which have larger er- 
rors), the second eigenmode (short-dashed) measures the dif- 
ference between collapsed and perpendicular configurations, 
this is sensitive to the dependence of Q3 on triangle shape. 




50 100 150 

a 



Figure 10. The three dominant, in terms of signal to noise [see 
Eq. (22] Q-eigenvectors for triangles with ri3 = 2ri2 = 24 Mpc/h 
from 100 HV mock catalogs. The signal-to-noise values are 9.1 
(solid), 8.4 (short-dashed), and 8.0 (long-dashed). 

the third eigenvector (long-dashed) gives the difference be- 
tween the configuration dependence at small and large a 
(down-weighted on the largest scales because of the large 
errors), this is sensitive to the scale dependence of Q3. The 
behavior of the eigenvectors shown in Fig. 10 is completely 
analogous to what was found in the bispectrum case (see 
Fig. 16 in Scoccimarro 2000). The condition in Eq. (20) is 
satisfied by 6 Q-eigenmodes (individual signal-to-noise val- 
ues are S/N = 9.1,8.4,8.0,6.7,4.9 and 1.1) in this particu- 
lar case with a total signal-to-noise of S/N = 17.0. Later, in 
section 5.5 and Figs. 12-14, we will illustrate the use of the 
Q-eigenmodes to recover the bias parameters. 

5.4 Gaussianity of the Likelihood Function 

A central assumption we made above is that the Q- 
eigenmodes are Gaussian distributed. This can and should 
always be checked against numerical realizations of the par- 
ticular survey geometry and scales being studied, otherwise 
the outcome of the likelihood analysis could be significantly 
biased, even if the covariance between different triangles is 
taken into account. Here we explore the validity of the Gaus- 
sian likelihood approximation for the three-point function, 
extending similar studies done for the bispectrum (Scoc- 
cimarro 2000) and higher-order moments (Szapudi et al. 
2000). 

As the survey volume increases, for a triangle of given 
scale and shape, the Gaussian approximation should become 
better as more such triangles contribute to its Q3 and al- 
though one is adding correlated quantities, one naively ex- 
pects to reach the central limit theorem if the correlations 
are small. However, it is not obvious this is enough be- 
cause at the same time the error bars for the bias parame- 
ters shrink, and thus a small systematic bias (if comparable 
to the small error bars) can still affect the determination 
of parameters. A simple rule of thumb is to compute the 
skewness of the distribution of Q3 and check whether this is 
much smaller than the shift in Q3 that results from displac- 
ing 61 , 62 by the error bars, if so the systematic error that 
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Figure 11. Histograms comparing the measured distribution of 3-point amplitudes Q3 (normalized to have unit variance) to the expected 
Gaussian distribution. The left panel corresponds to the 100 HV mocks, which results in a slightly skewed non-Gaussian distribution. 
The middle panel shows the distribution for the 5 dominant Q-eigenmodes Q using the same data, which is much more Gaussian. The 
right panel shows the first 5 dominant Q-eigenmodes Q in the 22 2dFGRS mocks. All cases correspond to triangles with fixed sides 
ri2 = 2ri3 = 12 Mpc//i and different values of a. 



results from assuming Gaussianity can perhaps be safely ig- 
nored (see discussion in Hui & Gaztaiiaga 1999, Scoccimarro 
2000 and Szapudi et al. 2000 for more details). 

Here we compare the histograms of values of Q3 in 
the galaxy mocks to a normalized Gaussian distribution: for 
each triangular shape, a, ri2 and ria, we estimate the mean 
value of Q3, Qa, and its variance ag in the mock catalogs. 
We then build a histogram of (Q3 — Qz)/(jq values for dif- 
ferent realiazations and also different configuration angles 
Q, but with fixed r\2 and ria sides. Although in general the 
distribution is not far from a Gaussian, it sometimes shows 
significant deviations. 

The left panel in Fig. 11 shows the worst situation we 
have encountered in our study. The distribution is somewhat 
skewed, although this is not a large effect (notice the loga- 
rthmic scale) . This is much more reassuring than the analo- 
gous result for IRAS surveys (see e.g. Fig. 15 in Scoccimarro 
2000) where the distribution becomes strongly non-Gaussian 
for smaller galaxy samples. In our case, the 2dFRGS or HV 
mocks are large enough to avoid these severe biases and the 
Gaussian likelihood seems a good approximation, but the 
distribution still has some small degree of skewness. In the 
left panel in Fig. 11 the dimensionless skewness of the dis- 
tribution is -1-1.1, which translates in a bias of about 10% in 
the estimated Q3. 

The relevant distribution in our method is not so 
much the distribution of Q3, but the distribution of the Q- 
eigenmodes. These are shown in the middle (HV mocks) and 
right (2dFGRS mocks) panels of Fig. 11 (notice the change 
in vertical scale from the left panel). The distribution seems 
better approximated by a Gaussian, which is expected since 
these are the highest signal to noise modes and should ap- 
proach Gaussianity fastest. The results correspond to trian- 
gles with fixed sides ri2 = 2ri3 = 12 Mpc//i, but similar 
results are found for other scales. The dimensionless skew- 
ness in both cases is of order 0.1, much smaller than in the 
direct distribution of Q3. This translates into 1% bias in Qa, 
which is smaller that the sampling variance from the single 
HOD mock. 

We conclude that a Gaussian likelihood provides a good 
approximation for our particular situation, specially when 
we express our results in terms of Q-eigenmodes. We will 



therefore use the Gaussian test in Eq. (16), restricted to 
the subspace of the dominant Q-eigenmodes from the SVD 
criterion, Eq. (20) 



5.5 Recovering Bias Parameters 

Figure 12 shows how well we can recover the values of 61 
and C2 using the measurements of Qa for ria = 2ri2 = 6, 12 
and 24/i~^Mpc assuming that the clustering pattern follows 
that of the Mr < —20 HOD mock galaxies. We use two 
different sets of mock catalogs (HV and 2dFGRS) to esti- 
mate the error bars and Q-eigenmodes to assess the sen- 
sitivity of each survey geometry. The mean values for the 
dark matter and galaxy Q3 are from the same simulation 
(VLS). From inside out, the values of the contours denote 
Ax^ = X^~Xmin = 2.4, 6.2, 11.8 away from the minimum 
value in Eq. (16), corresponding to 1, 2 and 3cr confidence 
intervals for a Ax^ distribution with two free parameters. 
Each panel also displays the total S/N, according to Eq.[22], 
added in quadrature for the dominant Q-eigenmodes used in 
each case. 

As can be seen in Fig. 12, there is a main direction 
that is best determined corresponding to the first eigen- 
mode, this is basically the direction perpendicular to the line 
(Qavg + C2)/6i = const, where Qavg is the average value of 
Qa. As small scales are approached, Qa develops the charac- 
teristic U-shape anisotropy and if the errors do not decrease 
fast enough (top panel) a near degeneracy develops in the 
b\ — C2 plane. The reason is that hi scales the shape depen- 
dence of (5a(a), while C2 only shifts the curve up or down. In 
the limit where the "U" is made of vertical and horizontal 
lines alone, scaling (&i) and shifting (02) produce equivalent 
(degenerate) effects. For large scales, the contours become 
wide because of the lack of enough independent triangles. As 
a result of these two limits, the best configurations for con- 
straining the b\ — C2 parameter space in 2dFGRS are those 
at intermediate scales (ria ~ 2ri2 — 12/i~^Mpc). 

The crossed circle in Fig. 12 corresponds to the fei, C2 ex- 
pected for the input galaxy mock catalogues using the HOD 
predictions in Eq. (14). These predictions are also shown as 
dashed lines in Fig. 8 and provide a very good description 
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Figure 12. Constraints on bi and C2 for triangles of different scales, as labeled in each panel. In all cases we compare measurements 
of Q-j in dark matter simulations to galaxy HOD mock catalogs. The top panels use errors (and Q-eigenmodes) from 2dFGRS mocks, 
while the bottom panels use HV mocks (smaller contours because HV have 2-3 times the effective volume of 2dFGRS mocks). 



of the mean bias relation, at least on large scales. Note that 
there are no free parameters here, we are not fitting for bi 
and C2, but rather checking we recover the values used for 
6i and C2 to create the mock catalogs using parameters in 
Table 3. Also recall that we are using a single HOD mock 
(from VLS simulations) to test the predictions, which means 
that there is some variance in the estimation and we should 
not expect the crossed circle on top of the best fit value. 

On smaller scales, our treatment of biasing breaks down 
because one can no longer use the expansion in Eq. (11) nor 
the simple HOD predictions based on Eq. (14). This failure 
can be seen in the bottom left panel of Fig. 12, where one 
finds significant (> 3cr) deviations from the predicted values. 

Figure 13 illustrates the importance of using the dom- 
inant Q-eigenmodes, specially when we only have a small 
number of mock catalogs or samples to estimate the co- 
variance matrix. As illustrated by the top right panel, if 
a noisy covariance matrix is inverted directly one recovers 
the incorrect values with unrealistically small error bars due 
to numerical instabilities. When we restrict the number of 
Q-eigenmodes the errors become larger (compare top and 
bottom left panels), but the estimation is not biased by the 
number of mocks used (compare two bottom panels). Note 
how the S/N is wrongly estimated in the top right panel 
and corresponds to a lower estimate when we use only the 
dominant Q-eigenmodes . 

Finally, Fig. 14 shows the constraints for different vol- 
ume limited samples for the 2dFGRS mocks. The case 



—21 < Mt < —20 provides the best compromise between 
sample volume and number density. The brighter sample has 
more volume, but is dominated by shot-noise. The fainter 
samples with higher density do not cover a large enough 
volume. These results are qualitatively similar to Fig. 9 in 
Sefusatti & Scoccimarro (2004). 



6 CONCLUSIONS 

We have studied the three-point function amplitude Q3 in 
dark matter and mock galaxies with particular emphasis on 
the effects of redshift-space distortions and galaxy biasing. 

Whereas the dependence of Q3 on triangle shape is 
partialy erased in redshift space on large scales ri2 > 
10/i~^Mpc, at smaller scales Q3 develops a characteristic 
U-shape anisotropy due to the effects of velocity dispersion. 
This is expected based on visual inspection of the distri- 
bution in redshift space (the finger of God effects). But 
nevertheless, it has not been detected before in simulations 
nor observations. We studied why this is so and concluded 
that low resolution in the configuration angle a can suppress 
this sharp feature due to mixing of collapsed configurations 
(which are few) with the more numerous (with much less 
amplitude) configurations away from a — 0, 180 deg. This 
problem can be exacerbated by using the v parameter as 
a probe of angle (instead of a), which has been a common 
choice in small-scale studies. Moreover, errors are largest for 
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Figure 13. Constraints on 61 and C2 using the data in Fig. 8 
(Mr < -20) at ri3 = 2ri2 = 24/x~lMpc. The input values (in 
Table 3) are denoted with a crossed circle. The top panels use 
all the Q-eigenmodes (equivalent to a direct inversion of the co- 
variance matrix), while the bottom panels use only the dominant 
Q-eigenmodes (Eq. [20]). Left panels use 100 HV mocks while right 
panels use only the first 22 HV mocks. Using the direct covariance 
inversion with too few mock catalogs (top right panel) typically 
gives systematically biased results and unrealistically small error- 
bars and large S/N. 
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Figure 14. Constraints on 61 and C2 from realistic 2dFGRS 
mocks (using the dominant Q-eigenmodes) of different absolute 
magnitudes (the volume-limited samples in Table 2). Here we only 
use triangles with ri3 = 2ri2 = 12/i~-'^Mpc. 



archical ansatz at small scales (less than 6 Mpc/h), much 
less so than in real-space, contrary to statements in the lit- 
erature. On the other hand, the scale dependence of Q3 is 
weakly suppressed in redshift-space when compared to real- 
space, thus if the U-shape anisotropy of Qa is missed due to 
poor resolution, a hierarchical three-point function becomes 
a good model. 

We presented a detailed method for obtaining con- 
straints on galaxy bias parameters from measurements of 
Qz in current galaxy redshift surveys, based on the Q- 
eigenmode analysis developed for the bispectrum (Scocci- 
marro 2000). We showed that the covariance matrix of Q3 
has non negligible extra diagonal components and presented 
an SVD inversion method that can be stable even if too few 
mock catalogs are used for the computation of the covari- 
ance matrix, with the only drawback in this case being an 
overestimate of the error bars. We studied the Gaussianity of 
the likelihood function and showed that for surveys such as 
2dfGRS Gaussianity is a good approximation for the distri- 
bution of the Q-eigenmodes, although some small skewness 
is present. 

We discussed the sensitivity of Qz to bias parameters 
as a function of the scale included in the analysis, and pre- 
sented results for how well we can expect to recover the lin- 
ear and quadratic bias parameters in the case of 2dFGRS. 
We showed that our method recovers the bias parameters 
introduced into mock galaxies by an HOD prescription and 
is also able to handle potential systematics in case a smaller 
number than ideal of mock catalogs is used to estimate the 
covariance matrix. 
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the collapsed configurations, which both makes it more dif- 
ficult to detect and is subject to larger estimation biases. 

Our detection of the U-shape anisotropy is completely 
analogous to the results from the small-scale bispectrum, 
which shows a very similar effect (Scoccimarro et al. 1999), 
and we see it in every cosmological model and including 
galaxy bias (see Figs. 3, 8). We therefore conclude that the 
redshift-space three-point function does not follow the hier- 
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